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/-s ABSTRACT 

CNJ Recent observation of the tidally-excited stellar oscillations in the main-sequence 

binary KOI-54 by the KEPLER satellite provides a unique opportunity for studying 
dynamical tides in eccentric binary systems. We develop a general theory of tidal 
excitation of oscillation modes of rotating binary stars, and apply our theory to tidally 

^^ excited gravity modes (g-modes) in KOI-54. The strongest observed oscillations, which 

C* . occur at 90 and 91 times the orbital frequency, are likely due to prograde m — 2 modes 

^ I (relative to the stellar spin axis) locked in resonance with the orbit. The remaining 

flux oscillations with frequencies that are integer multiples of the orbital frequency are 
likely due to nearly resonant m = g-modes; such axisymmetric modes generate larger 
flux variations compared to the m = 2 modes, assuming that the spin inclination angle 
jrt of the star is comparable to the orbital inclination angle. We examine the process of 

^^ resonance mode locking under the combined effects of dynamical tides on the stellar 

spin and orbit and the intrinsic stellar spindown. We show that KOI-54 can naturally 
evolve into a state in which at least one to = 2 mode is locked in resonance with 
the orbital frequency. Our analysis provides an explanation for the fact that only 
oscillations with frequencies less than 90-100 times the orbital frequency arc observed. 
We have also found evidence from the published KEPLER result that three-mode 

■^+ nonlinear coupling occurs in the KOI-54 system. We suggest that such nonlinear mode 

coupling may explain the observed oscillations that are not harmonics of the orbital 
frequency. 
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d 1 INTRODUCTION al. (2003), Ivanov & Papaloizou (2004) and Papaloizou & 

^. , , . , . , . , . Ivanov (2010). 
iides play an important role m binary star systems and m 

star-planet systems. While numerous studies of tidal effects Recent observations of the binary star system HD 

have been based on the so-called equilibrium tide theory, 187091 (KOI-54) by the ifepfer satellite provide a unique op- 

which parametrizes tidal dissipation by an effective tidal lag portunity for studying dynamical tides in eccentric binaries. 

angle/time or tidal quality factor (e.g., Darwin 1880; Goldre- KOI-54 consists of two A stars (mass Mi,2 = 2.32, 2.38 M© 

ich&Soter 1966; Alexander 1973; Hut 1981), the underlying and radius i?i,2 = 2.19, 2.33 i?©) in an eccentric (e = 

physics of tides in fluid stars and planets involves dynamical 0.8342) orbit with period P = 41.805 days (Welsh et 

excitations of waves and oscillations by the tidal force (see al. 2011). The binary is nearly face-on with orbital inclina- 

Ogilvie & Lin 2007 and Zahn 2008 for recent reviews). Tides tion iorb = 5.52°. In addition to periodic brightening events 

in highly eccentric systems are particularly rich in their dy- caused by tidal distortion and irradiation of the two stars 

namical behavior, since wave modes with a wide range of during their close periastron passages, the power spectrum of 

frequencies can be excited and participate in the tidal in- the Kepler light curve revealed 30 significant (with a signal- 

teraction. Various aspects of dynamical tides in eccentric to-noise ratio ^ 7) stellar pulsation modes. The observed 

binaries have been studied by Lai (1996,1997), Kumar & mode periods range from 45 hours to 11 hours, correspond- 

Quataert (1998), Witte & Savonije (1999,2001), Willems et ing to the mode frequency fa ranging from 22.42 /orb to 

91 /orb (where /orb ~ P^^ is the orbital frequency). Most 

interestingly, twenty-one of these mode frequencies are inte- 

* Email: derg@astro.cornell.edu; dong@astro.cornell.edu ger multiples of /orb (with the ratio fa/ foih differing from an 
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integer by 0.01 or less). The two dominant modes have fre- 
quencies that are exactly 90 and 91 times /orb, with the cor- 
responding flux amplitudes of 297.7 /imag and 229.4 /imag, 
respectively. 

While dynamical tides in massive-star binaries have 
been studied before (e.g., Zahn 1977, Goldreich & Nicholson 
1989 for circular binaries; Lai 1996,1997, Kumar & Quataert 
1998 and Witte & Savonije 1999,2011 for eccentric binaries), 
KOI-54 represents the first example where tidally excited 
oscillations are directly observed and therefore serves as an 
explict demonstration of dynamical tides at work in the sys- 
tem. As discussed in Welsh et al. (2011), the observed os- 
cillation modes are puzzling: over 20 of the observed modes 
are nearly exact integer multiples of the orbital frequency, 
yet several others are not. It is not clear why the dominant 
modes are so prominent, e.g., why modes with frequencies 
90, 91, 44, and 40 times /orb are clearly visible, and yet 
modes with frequencies greater than 91/orb and those less 
than 20/orb appear to be absent. 

The goal of this paper is to explain some of the ob- 
servational puzzles related to KOI-54 and to develop the 
general theoretical framework for studying tidally-excited 
oscillations in eccentric binary systems. 

Our paper is organized as follows. In Section 2, we derive 
the general equations for calculating the energies of tidally 
excited oscillation modes in an eccentric binary. Our theory 
improves upon previous (and less rigorous) works, and pro- 
vides a clear relationship between the resonant mode energy 
and non-resonant mode energy. In Section 3 we study the 
properties of non-radial g-modes relevant to the stars in the 
KOI-54 system and calculate the non-resonant mode ener- 
gies - these serve as a benchmark for examining the effect 
of resonances. In Section 4 we present our calculations of 
the flux variation due to tidally-forced oscillations. We show 
that the observed flux variation in KOI-54 can be largely 
explained when a high-frequency mode is locked into reso- 
nance (with the mode frequency equal to 90/orb). In Section 
5 we study the possibility of resonance locking. We show 
that the combination of the secular tidal orbital/spin evo- 
lution and the intrinsic spindown of the star (e.g., due to 
stellar evolution) may naturally lead to resonance locking of 
a particular mode. Our analysis demonstrates that (in the 
KOI-54 system) a mode with frequency around 90/orb can 
be resonantly locked, while modes with higher frequencies 
cannot. In Section 6 we discuss the origin of the observed 
modes in KOI-54 with frequencies that are not an integer 
multiple of /orb, including the evidence of nonlinear mode 
coupling. We conclude in Section 7 with a discussion of fu- 
ture prospects and remaining puzzles. 



2 DYNAMICAL TIDE IN ECCENTRIC 
BINARY STARS: GENERAL THEORY 

We consider the tidally-excited oscillations of the primary 
star of mass M and radius R by the companion of mass M'. 
The gravitational potential produced by M' can be written 



where r^ = (r, 9,<j)i — (p + Q,st) is the position vector (in 
spherical coordinates) relative to the center of star Af (the 
azimuthal angle (p is measured in the rotating frame of the 
star, with the rotation rate fls and the spin axis aligned 
with the orbital angular momentum), D{t) is the binary 
separation and $ is the orbital true anomaly. The dominant 
terms have I = |77i| — 2 and / = 2, m = 0, and for these 
terms W2±2 = (STr/lO)^/^ and W20 = (7r/5)^/^ The linear 
response of star M is specified by the Lagrangian displace- 
ment ^(r,i), which satisfies the equation of motion (in the 
rotating frame of the star) 






+ 2a,s 



di 



x- + C.^ = -Vt/, 



(2) 



where C is a self-adjoint operator (a function of the pressure 
and gravity perturbations) acting on ^ (see, e.g., Friedman 
& Schutz 1978). A free mode of frequency Wq (in the rotating 



frame) with ^„(r,t) = ^^(r) e~ 



im4> — iijJcit 



satisfies 



- ujU^ - 2iujMs X ^„ -h C • ^„ = 0, 



(3) 



where {a} denotes the mode index, which includes the radial 
mode number n, the polar mode number L (which reduces 
to I for spherical stars) and the azimuthal mode number 
m. We carry out phase space mode expansion (Schenk et 
al. 2002) 



d£,/dt 



Y,c^{t) 






(4) 



where the sum includes not only mode indices, but also both 
positive and negative uja ■ Note that the usual mode decom- 
position, ^ = "^^a' '^'^^a (with the sum including only mode 
indices), adopted in many previous studies (e.g., Lai 1997, 
Kumar & Quataert 1998; Witte & Savonije 1999), are rig- 
orously valid only for non-rotating stars. Using the orthog- 
onality relation {£,^,2ifl^ x C'> + i^c, + ^a'ii^c^^a') = 
(for a / a'V where {A, B) = Jd^xp{A* ■ B), we find (Lai 
& Wu 2006^] 

C„ + {iuJcc +'^cc)Ca = - — (^„(r),-VC/) 

where 7a is the mode (amplitude) damping rate, and 

Q^ = {i^,\7{r'Yi^)), (6) 

ec,=wc + {i^,in,x^^), (7) 

and we have used the normalization {^cuia) ~ 1- The quan- 
tity Qa (called the "tidal overlap integral" or "tidal coupling 
coefficient" ) directly relates to the tidally excited mode am- 
plitude. We shall focus on I = 2, m — and |m| = 2 modes 
in the following (although we will continue to use the no- 
tations I and m so that it would be easy to generalize to 
high-order tides). 

The general solution equation ([5| is 



Ca{t) 



iGM'Wi„,Qc 
2eo.Di+^ ' 



i<ya^-\-"'la t — im^ 



dt, (8) 



C/(r.,t) = -GM'J2^^^ e-™*'*) 



^ As noted before, in this paper we restrict to aligned spin-orbit 
Yi (0 (b) (1) configurations for simplicity. Generalization to misaligned sys- 

tems is straightforward (Lai & Wu 2006; see also Ho & Lai 1999). 
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assuming Ca{to) — 0, where 



: iUa + mQs 



(9) 



is the mode frequency in the inertial frame. Let tj — (2j — 
l)P/2 (with J = 0, 1, 2, • • • ) be the times at apastron. After 
the fcth periastron passage, the mode amplitude becomes 



Co.{tk) =(Ac,)e''""=*'=-('"°+-^°>^/' 



-(icra+'-fa)kP 



-(i<ya+la)P 



with 



Acc = 



P/2 



iGM'Wi^Q 



dt „, , " — e 

P/2 2eai?'+i 



oi icTct +'ya£ — imcl" 



(10) 



(11) 



For 'jakP S> 1, the steady-state mode energy in the inertial 
frame 



±ya — ^^a&a\(^ci 



(12) 



becomes (Lai 1997) Q 



A£;„ 



2 (cosh 7a P — cosac-P) 

AS, 



(13) 



4sin2(a,P/2) + (7„P)2' 

where the second equality assumes 7a P <^ 1. Here A^Eq 
is the energy transfer to the mode in the "first" periastron 
passage: 

AiJ„ = ^^(^l ^^^IQ.i^i™(aJ^ (14) 



R 



where Dp = a{\ — e) is the periastron distance (a is the 
orbital semi-major axis) and 



KlmiOa) — 



Wvrn 

27r 



p/2 



dt 



p/2 



i) 



l+i 



e ° . (15) 



Note that in equation ( |14[ |, both Q^ and 7^;^ are dimen- 
sionless (in units such that G = M = R — 1). 

Equation ( 13 1 shows that when aaP is not close to 2nN 



(where A*' is an integer), the steady-state mode energy is ap- 
proximately AEa ■ Thus AEa serves as a benchmark for the 



non-resonant mode energy. Equation ( 13 1 provides a simple 



relationship between the actual mode energy Ea and the 
non-resonant mode energy AEa- 



3 STELLAR OSCILLATION MODES AND 
NON-RESONANT MODE ENERGIES 

We construct an M = 2.35 Mq stellar model using the 
MESA code (Paxton et al. 2010). We assume solar metallic- 
ity and evolve the star until its radius reaches R — 2.34 Pq. 
These parameters are close to star A/2 in KOI-54. Figure [l] 
displays a propagation diagram for our stellar model. The 
star has a small convective core inside radius r — 0.09P. We 
make sure that the stellar model has thermodynamically 
consistent pressure, density, sound speed and Brunt- Vasala 
frequency profiles. We have computed the I = 2 adiabatic 

■^ Equation (1131 was derived in Lai (1997) in an approximate 
manner (since mode decomposition was not done rigorously), and 
physical arguments were used to get rid of a fictitious term. 



g-modes for this non-rotating stellar model, including uia, 
Qa, and the mode mass Ma = (^^ ' ^a)/l^a(-R)l^- Here, the 
magnitude of the surface displacement is defined by 



\UR)\ 



dn^^iR)-CiR) 



= &l{R) + Ki + m±liR), 



(16) 



where the r and _L subscripts denote the radial and horizon- 
tal components of the displacement vector, respectively. 

We use equation (141 to compute the non-resonant 
mode energy E^ ~ AEa- The corresponding surface dis- 
placement £,a{R) is then obtained from 



\iaiR)\ 



Ea 



1/2 



(17) 



Figure [2] shows the energy AEa and surface amplitude 
of the radial component of the displacement, £,ra{R)j of 
tidally excited modes away from resonances for the KOI-54 
parameters. The radial displacement is directly related to 
the flux variation due to the oscillation mode (see Section 
4). The most energetic modes have frequencies 80f2 ^ (Tq, <J 
14057, where Q is the orbital angular frequency, depending 
on the stellar spin rate and the value of m. Low-order (high 
frequency) modes have larger values of Qa but have smaller 
values of Kim{<^a), so medium-order (n ~ 15) modes have 
the largest values of AEa rl 

Figure [2] also shows the magnitude of the displacement 
|4^(P)| for m = 2 modes in the zero spin limit. The total 
displacement is much larger than the radial displacement for 
low-frequency modes because these modes are characterized 
by large horizontal displacements and are concentrated near 
the surface. Consequently, these modes have lower mode 
mass and the maximum of |^q(P)! shifts to lower frequencies 
compared to the maximum of AEa ■ 

The rotation rates of the KOI-54 stars are unknown. 
Spectroscopic observations constrain VJ-ot sinig ^ 10 km s^^, 
corresponding to K-ot ^ 100 km s~^ if the spin inclination 
angle is equals is = 5.5° (Welsh et al. 2011). This implies 
that the spin period Ps <: 1.2 days and that Qs SS 30Q,. 
Although the classical equilibrium tide theory (e.g.. Hut 
1981) is not expected to be valid for our system (Lai 1997), 
we can adopt the pseudosynchronous rotation frequency 
f7s = fips = 16.5f2 [corresponding to fps — 1/(2.53 days)] as 
a fiducial value. 

To account for the effect of stellar rotation on the tidally 
excited modes, we adopt the perturbative approximation, 
valid when Qs is less than wi (the mode frequency in the 
zero-rotation limit). The mode wave functions and Qa are 
unchanged by the stellar rotation, while the mode frequen- 
cies are modified according to cj^ = uja — mCniQ,s and 
aa = (jJa + m{l — C„i)Qs, where Cni > is a constant (e.g. 
Unno et al. 1989) — our calculation gives Cni ~ 0.16 for all 
relevant modes. Note that in this approximation, Sa = ^a ■ 
More accurate results can be obtained using the method of 
Lai (1997). 

Stellar rotation increases the inertial-frame frequency 



^ The orbital frequency at periastron is fp = /orb(l + e)^''^/(l " 
e)3/2 = 20.06 /orb = 1/(2.084 d), thus m = 2 modes with era /r2 ^ 



40 have the largest values of Ki^. For m ■■ 
1 have the largest values of Ki^. 



■ 0, modes with aa/^ > 
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Figure 1. Propagation diagram for our M = 2.35A/0, R = 2.34/?0 stellar model, showing the value of A' (solid black line) and the 
Lamb frequency, L2 (dashed black line), in units of (GM/R^)^'^. The horizontal red lines denote important angular frequencies for the 
KOI-54 system, including (from top) the dynamical frequency of the star, (GM/iJ^)^ "; the highest frequency mode observed in KOI-54 
((To, = 9in); the lowest frequency mode observed (ctq = 22.42f2); and the orbital angular frequency, Q. The y— axis on the right-hand 
side displays the corresponding periods, in units of hours. 



act of the m, > modes, causing the dominant modes to have 
higher observed frequencies. However, rotation also shifts 
the tidal response to higher order g-modes with smaller val- 
ues of Qa , so we expect rotation to lower the mode energies. 
Figure [2] confirms that finite (prograde) stellar rotation in- 
deed shifts the dominant mode energy and surface ampli- 
tude to higher-order g-modes, giving rise to smaller AEa 
and ^ra{R)- 

The somewhat erratic features of AEa and ^q as seen 
in Figure [2] for high-order (low frequency) modes are due 
to mode trapping effects created by the thin sub-surface 
convection zone in our stellar model. For these high-order 
modes, the tidal overlap integrals (Qa) depend on the pre- 
cise shape of the mode wave function, and they can be easily 
affected by the detailed properties of the stellar envelope. 
Care must be taken in order to obtain reliable tidal overlap 
integrals for these high-order modes (see Fuller & Lai 2010). 

The damping rate of a mode, 7^, can be estimated in 
the quasi-adiabatic limit via 



the local radial wave number, and r+ and r_ are the bound- 



J r_ 



(18) 



aries of the mode's propagation cavity. Equation ( 18 1 can 



be derived in the WKB limit from the quasi-adiabatic work 
function of a mode (Unno et al. 1989).rjThe bottom panel of 
Figure [2] shows the values of ja calculated via equation ( 18 1 



for the modes of our stellar model. Lower frequency g-modes 
have larger damping rates due to their larger wavenumbers 
and because they propagate closer to the surface of the star 
where the diffusivity is larger. Equation ( 18 1 provides an 



estimate for mode damping rates via radiative diffusion in 
the quasi-adiabatic limit; however, fully adiabatic oscilla- 
tion equations must be solved for low-frequency modes for 
which 7a, is comparable to the mode frequency. Furthermore, 
modes may also damp via nonlinear processes (see Section 



where x is the thermal diffusivity, k^ ~ l{l + 1)N'^ /{r^uP') is 



■* While our paper was under review, we became aware of the 
paper by Burkart et al. (2011), who also used equation l |18| to 
estimate the mode damping rate. The bottom panel of Figure [2] 
was added after we saw the Burkart et al. paper. 
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Figure 2. Non-resonant mode energy AEa (top), amplitude of the radial surface displacement Ci-a(^) (middle), and mode damping 
rates -faP (bottom) as a function of mode frequency era (in units of the orbital frequency Q). The I = m = 2 modes are plotted as open 
circles for fig = and filled circles for fig = 16. 5f!, while m = modes are plotted as triangles. The middle panel also displays the 
magnitude of the surface displacement, |^a(fl)|, for the I = m = 2 modes with fis = (open circles connected by lines). The highest 
frequency modes shown are the gg mode for the m = and non-spinning m = 2 cases, and the gi2 mode for the spinning m = 2 case. 
The mode energy is in units of GM^ /R and the displacement is in units of R. 



4 FLUX VARIATION DUE TO 

TIDALLY-FORCED OSCILLATIONS 

In section [2l we considered the tidal response to be com- 
posed of the sum of the star's natural oscillation modes, 
each having its own steady-state energy Ea. This provides 
a simple relation between the resonant mode energy and 
non-resonant energy. The observed magnitude variation of 
KOI-54 has over 20 components with frequencies that are al- 
most exact multiples of the orbital frequency, i.e., they have 
o"a = (A'^ + e)f2 with |e| ^ 0.01. In this section, we examine 



the tidal response as a sum of oscillations at exact multiples 
of the orbital frequency. Each component Uim of the tidal 
potential can be decomposed as 



Ulm= ^n Ylm[0,(pt) 2_^ FNme , (19) 

N—~oo 



where FNm. is defined by the expansion 



,dJ ^" 



J2 p«^ 



(20) 
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and is given by 

-FjVm ~ 



d* 



COS [iV(* - esin*) - m<^(t)\ 
(l - ecos*)' 



(21) 



Here, a is the semi-major axis of the orbit, and '^ is the 
eccentric anomaly. Note that Fjvm is related to Kim [see 
equation ^\ by iC^^K = NO.) = Wim{l-ef+^Q.-^FNm. 
Inserting equation ( 19 \ into equation ([5| yields 



+ (iuJa + ')a)c 



iGM'Wlr, 



2ec.a'+^ 



y^ -FiVme 



i{mns—N(l)t 



whose non-homogeneous solution is 

GM'WimQc ^ FNme 



«(i) 



2e„a'+i 



E 



(jvn-mns)t 

((Ta - iVri) - i'Ja 



(22) 



(23) 



«.(r) 
terms of r^ 
frame, with 



Amcp 



The total tidal response is ^(r,i) = 51^^ Ca(t)^Q,(r), where 
When the displacements are expressed in 
{r,6,(j}i) (the position vector in the inertial 
= (f) + mQgt), we have 






N — ~-oo a 



^ ^ GM'W,„0, 



Z.E 

FNme~^ 



2£ca'+i 









F- 



iVme 



(ctq - NO.) - i7„ (ctq + NO.) - i"/a 



+ c.c. (24) 



where c.c. denotes the complex conjugate, and "^^a' unplies 
that the sum is restricted to modes with cjq, > (includ- 
ing both m > and m < modes) (by contrast, "^^ in- 
cludes both (jja > and uja < terms, as well as both 
m > and m < 0). We have omitted the A^ = term for 
simplicity because it is not part of the dynamical response. 
Each ^jy (oscillating at frequency NQ) is then a sum over 
the star's oscillation modes, with large contributions coming 
from nearly-resonant modes with aa « ATi. 

We use the result of Buta & Smith (1979) (see also 
Robinson et al. 1982) to estimate the magnitude variation 
(Amag) associated with each component of the tidal re- 
sponse. The magnitude variation of the star has two pri- 
mary components: a geometrical component, (Amag)^, due 
to distortions in the shape of the star, and a temperature 
component, (Amag)T, due to perturbations in the surface 
temperature of the star's photosphere. 

For an / = m, = 2 mode with surface radial displace- 
ment (,ra{R), the amplitudes of the bolometric magnitude 
variations are 

(m=2) . . ^ „ _ „.__2 _■ r2 



(Amag)^"^ 



-1.77isin is 



l{l- 



Ta 



R 



and 



/A N (tti — 2) 

(Amag)^e 



-0.42(a,-f/30sin^, ^^. 



(25) 



(26) 



Here, T2 — 5/3 is the adiabatic index of gas at the surface 
of the star, Cj = uja/ \fGMjR? is the dimensionless mode 



frequency, and 7; « 0.3 and a; + /3; ~ — 1.2 are bolomet- 
ric limb darkening coefficients appropriate for A stars. For 
KOI-54, the value of the spin inclination angle is (the angle 
between the spin axis and the line of sight) is unknown, but 
we may use the system's orbital inclination of Jorb = 5.5° as 
a first guess, although it is possible that the star's spin axis 
is inclined relative to the orbital angular momentum axis. 
Similarly, for / = 2, m = modes, the amplitudes of the 
magnitude variations are 



(Amag)[j5' "' — -1-47! (3cos^i3 



l{l + l) 



1 



V2 



r2 

^r^{R) 



R 



and 



(Amag) 



(m — 
aG 



-0.34 (a, + A) (3 ( 



1)' 



R ■ 



(27) 



(28) 



In equations (251 and (27 1, the factor of [/(/-|-l)/cij„— 4- 



Cjc,] arises from the outer boundary condition formulated by 
Baker & Kippenhahn (1965) (see also Dziembowski 1971), 
in which the radial derivative of the Lagrangian pressure 
perturbation vanishes at the outer boundary (in our case, 
the stellar photosphere) , such that the mode energy vanishes 
at infinity in an isothermal atmosphere. It is not clear how 
well this boundary condition applies for a real star, nor how 
it should be modified when non-adiabatic effects are taken 
into account. For the g-modes considered here, the factor 
l{l+l) / ujIi S> 1, causing the temperature effect to overwhelm 
the geometrical effect for magnitude variations. 

The magnitude variation for each frequency A''^ is a 
sum of the variations due to individual modes: 

GM'WlmQc.a^i,r^{R) 



(Amag)iv = ^ 



FnttiC 



+ 



2eca'+i 



-NmC 



{ffa — Nil) — i'ya ((7a + Nfl) — ija 



R 



+ C.C. 



(29) 



where a^ is the constant in front of (,ra{R)/R in equations 



(25l-(28l, andiV>0. 



Using the values of Sa, Qa, ^a(^)i ^^'^ 7q calculated 
in Section [2) w e compute each term (Amag)]v in equation 
( 29 1 . Figure |3| shows a plot of the magnitude variation as a 
function of TV, along with the observed magnitude variations 
in KOI-54 as determined by Welsh et al. (2011). To make 
this plot, we have subtracted out the contribution from the 
"equilibrium tide", because the equilibrium tide is respon- 
sible for the periodic brightening of KOI-54 near periastron 
and was subtracted out by Welsh et al. (2011) in order to 
obtain the magnitude variations due to tidally-induced stel- 
lar oscillations. The equilibrium tide can be computed by 
taking the limit aa ^ NQ, i.e., by setting NO, — for each 
term inside the sum in equation ( |29| |. We adopt a spin in- 
clination angle is = 10°. 

Figure [S] includes the magnitude variations due to the 
temperature effect and those due to geometrical effects. For 
the adopted spin inclination (10°), the m = modes domi- 
nate the observed magnitude variations, although a nearly- 
resonant m = 2 mode can be visible. However, it is obvi- 
ous from Figure |3] that our computed variations due to the 
temperature effect of m = modes are appreciably larger 
than those observed. In fact, the magnitude variations due to 
geometrical effects reproduce the observed variations much 
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Figure 3. Magnitude variation Amag^y as a function of N due to the temperature effect (top) and geometrical effect (bottom) for 
I = m = 2 modes (filled circles) and I = 2, m = modes (triangles). We also plot the observed magnitude variations that are an integer 
multiple of the orbital frequency (red plus symbols), and those that are not an integer multiple of the orbital frequency (red x symbols). 
The plot uses ig = 10° and Qg ~ 15^2. 



more accurately. The over-prediction of the magnitude vari- 
ations from the temperature effect is mostly likely due to 
non-adiabatic effects in the stellar atmosphere, which ren- 
ders equations ( |25[ ) and ( |27[ ) inaccurate for low-frequency 
modes. Indeed, Buta & Smith (1979) also found that for 
main sequence B stars, the predicted magnitude variations 
due to the temperature effect for low-frequency modes were 
much larger than what was observed, and they speculated 
that the mismatch was due to non-adiabatic effects in the 
outer layers of the star. 

Our results depicted in Figure [3] suggest that most of 
the observed magnitude variations (with the exception the 
highest-frequency modes) in KOI-54 are due the geometrical 
surface distortions produced by m = modes that happen to 
be nearly resonant with a harmonic of the orbital frequency. 



The actual magnitude variations due to the temperature 
effect (when non-adiabatic effects are properly taken into 
account near the stellar photosphere) should not be much 
larger than the geometric effect. Note that in Figure [3] we 
have plotted the absolute values of the magnitude variations 
for the temperature and geometrical effects, but that the 
temperature and geometrical effects have opposite sign and 
tend to cancel each other out. 

To produce the highest-frequency modes (90/orb and 
91/orb) observed in KIO-54, in Figure [slwe have fine-tuned 
the spin of the star to fis ~ 150 such that the A*' = 90 
oscillation is nearly resonant with one of the star's oscilla- 
tion modes. In section 5 we argue that the A'' = 90 and 
A*' = 91 oscillations are likely due to an m, = 2 mode in each 
star that is locked in resonance. The fine-tuning produces 
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a magnitude variation due to the temperature effect that 
is comparable to one of the observed oscillations, although 
the predicted geometrical magnitude variation is somewhat 
less than what is observed. If the A^ = 90 and A^ = 91 os- 
cillations are indeed due to resonant m = 2 modes, then 
it is likely that the predicted magnitude variation due to 
the temperature effect (based on adiabatic approximation) 
is more accurate for these modes due to their higher frequen- 
cies. This is not unreasonable since non-adiabatic effects in 
the stellar photosphere are expected to be less important 
for low-order (high frequency) modes, a point also empha- 
sized by Buta & Smith (1979). In any case, a more in-depth 
analysis of the non-adiabatic oscillation modes (particularly 
their flux variations) for the stars of the KOI-54 system is 
needed. 



5 SECULAR SPIN-ORBIT EVOLUTION AND 
RESONANCE LOCKING 

In the previous sections we proposed that a resonance with 
ac, ~ 90n creates the largest observed magnitude variations 
in KOI-54. Here we study the secular evolution of the stellar 
spin and binary orbit, and how resonances may naturally 
arise during the evolution. Several aspects of tidal resonance 
locking in massive-star binaries were previously explored by 
Witte & Savonije (1999, 2001). 

Dynamical tides lead to spin and orbital evolution, with 
the orbital energy and angular momentum evolving accord- 
ing to 



-Eorb ~ — 2_^ "ijaEa, Jorb = — 2_^ ^JaEa 



m 



(30) 



where we have used the fact that the mode angular momen- 
tum (in the inertial frame) is related to the mode energy by 
Ja = (m/aa)Ea. The orbital frequency fl therefore evolves 
as 






— , 



where 



U= [4sm\nac,/n) + {'y^Pf 



(31) 



(32) 



and t^^ specifies the orbital decay rate due to a single non- 
resonant mode a: 

l-Eorbl 



1271^70 f M' 
(l-e)6 \JY 



\QaKlr, 



(33) 



Using the KOI-54 parameters and assuming (Jq 
have 



1.9 X 10"^°7q 



10-4 



(34) 



Resonance occurs when a mode has frequency a^ — N^ for 
an integer value of N. 



5.1 Need for Resonance Locking 

We first consider how likely it is to observe a mode near 
resonance when a^ is held constant (i.e., it does not change 
during the period of resonance-crossing) . 

As can be seen from Figure [2] the m = 2 modes dom- 
inate tidal energy transfer in the KOI-54 system, yet (as 
evidenced from Figure |3|, most of the visible modes (except 
A*' = 90, 91) are m = modes. Consider a particular m = 
mode near resonance, with aa = {N + e)fl and |e| <^ 1. Since 
the mode contributes little to the tidal energy transfer, the 
probability of being close to resonance (|ej < eo) is simply 
Pros — 2eo. Figure [3] indicates that an tti = mode with 
25 Ji A'' <J 80 requires |eo| <; 0.1 to be visible [mode visibil- 
ity scales as |(Amag)j oc l/|e|]. In each star there are about 
forty m — modes in this frequency range, thus we should 
expect about 8 of these modes from each star to be visible, 
in rough agreement with Figure [S] and the observations. 

Now let us consider modes that significantly infiuence 
the tidal energy transfer (these include m = 2 modes such 
as the A' — 90, 91 modes, but may also include m = modes 
very near resonance - if they occur). For the KOI-54 system, 
the modes that contribute significantly to the tidal energy 
transfer have a^ <: 2f2p ~ 40f2, thus TV = aa/0. <; 40. Con- 
sider a particular mode near resonance, and suppose that 
the tidal energy transfer is dominated by the resonant mode 
(q). The orbital decay rate during the resonance is given 
by P/P = -t-^ [(Svre)^ -1- (7cP)^]"\ Thus the time that 
the system spends "in resonance" (|e| < eo) is (At)e(j = 
td«(2eo/A') [(27reo)V3+(7<.P)'] ~ ^■Khd^\tQ\'^ /{?,N). By 
contrast, the time it takes the orbit to evolve between res- 
onances (for the same mode) is (Ai)nr ~ td/N, where td is 
the orbital decay timescale due to all non-resonant modes {td 
may be a factor of ~ 10 smaller than tda)- The probability 
of seeing a mode very near resonance is thus 



(Ml 

{At), 



Btt^ tda , 
3 td 



m 



(35) 



Figure 13] indicates that (for in inclination of is ~ 10°) we 
require jeo! ^ 10~^ for an m — 2 mode to be visible, for 
which the probability is Pres ^ 3x lO"**. Thus, at first glance, 
the chances of observing even a single tidally-excited m — 
2 mode in the KOI-54 system are shm. The A'' = 90,91 
modes require |eo| ^ 10~^, even if they are produced by m = 
modes. It is therefore extremely unlikely to observe the 
system with such large amplitude modes, unless (by some 
mechanism) they are locked into resonance. In the sections 
that follow, we outline such a resonance locking mechanism 
and how it applies to the KOI-54 system. 



5.2 Critical Resonance-Locking Mode 

We now consider the possibility of resonance locking for the 
m = 2 modes. Tidal angular momentum transfer to the star 
increases the stellar spin Qa, thereby changing the mode fre- 
quency aa- There exists a particular resonance, aa = Nci^, 
for which a^ and SI change at the same rate, i.e.. 






(36) 



so that the resonance can be maintained once it is reached. 
With Js = —Jorb and the mean stellar rotation rate 



Dynamical Tides in KOI- 54 9 



^s ~ Js/I (where I is the moment of inertia), we have 
&ci = mBaJs/I, where B^ = {doa/dCls)/m. Assum- 
ing that a single resonant mode dominates the tidal en- 
ergy and angular momentum transfer, we find da /a a = 
-m^BaEorh/iIN'^0.'^)- Thus 



N ) [n 



where 



Nc = m 



BajJ-a^ 
3/ 



2\ 1/2 



(37) 



(38) 



(with fi — MM' /Mt the reduced mass of the binary) corre- 
sponds to the critical resonance-locking mode. The subscript 
"tide" in equation ( |37[ ) serves as a reminder that the changes 
of (Ta and fi are due to the tidal interaction. 

For the KOI-54 parameters (with a ~ 397? and m = 2), 
we have 



N^ ~ 131 



0.8' 



1/2 



Vo.o 



0.05 



-1/2 



(39) 



where k — I /{MB?). If we assume that the star maintains 
rigid-body rotation during tidal spin up, then Be, = 1 — 
Cni- Numerical calculation for our stellar model gives Ba — 
0.84 (with very weak dependence on modes) and k ~ 0.040, 
giving A^c — 146; a less evolved star would have k ~ 0.047, 
giving Nc — 135. Note that in the above we consider only 
tides in one star (M) - when tides in the other star are 
considered, f2 would be larger and the effective value of Nc 
would be reduced. If identical resonances occur in the two 
stars, the effective value of Nc would be reduced by a factor 
of yM to Nc ^ 92 (for k = 0.05). See section 5.3.3 and 
equations ( |70[ )-(|7l| below. 

The above consideration assumes that the spin evolu- 
tion of the star is entirely driven by tide. In reality, the 
star can experience intrinsic spindown, either due to a mag- 
netic wind or due to radius expansion associated with stellar 
evolution. We may write tig = (f2s)tidc — l^slsd, where the 
second term denotes the contribution due to the intrinsic 
spindown. The corresponding rate of change for Ua is then 



Oa 
Oc 



(40) 



On the other hand, the orbital evolution remains solely 
driven by tides: 






(41) 



Comparing equations ( |40[ ) and (WTl, we see that Nc repre- 
sents the upper boundary for resonance locking: for N > Nc, 
resonance locking {a a jo a ~ Q/Q) is not possible. For 
N < Nc, resonance locking can be achieved: the value of 
{a-c,/o-a)tida depcuds on the closeness to the resonance, which 
in turn depends on the intrinsic spindown timescale of the 
star [see equations ( |54[ ) and ( |55| l below]. 

The fact that Nc naturally falls in the range close to 
N = 90 and A'^ — 91, the most prominent modes observed 
in KOI-54, is encouraging. In the next subsection, we con- 
sider how the system may naturally evolve into a resonance- 
locking state with N < Nc- 



5.3 Evolution Toward Resonance 

As noted above, even in the absence of tides, the star can 
experience spin-down, either due to magnetic wind or due 
to radius expansion associated with stellar evolution. Fur- 
thermore, the intrinsic mode frequencies will change as the 
internal structure of the star changes due to stellar evolu- 
tion. We now show that these effects can naturally lead to 
evolution of the system toward resonance locking. 

The evolution equation for the stellar spin reads Cls — 
— {Joib/I) — {^s/tsd), where tsd is the "intrinsic" stellar spin- 
down time scale associated with radius expansion and/or 
magnetic braking. Using equation (|30[), we have 






= E 






n' 



n. 






1 

isd 



(42) 



From equation ( 30 1 , we find that the orbital frequency and 



eccentricity evolve according to 






E 



tdo 



IE 






mil 



cr„(l-e2)l/2 



(43) 



(44) 



To leading order in Qs, the mode frequency depends on f^s 
via 



,(0) 



4- mBafls 



(45) 



For simplicity, we assume that the secular change of aa is 
only caused by the fis-evolution (e.g., we neglect the vari- 
ation of uja due to stellar evolution - this effect can be 
absorbed into the spindown effect on the mode; see Witte & 
Savonije 1999). 



5.3.1 Single Mode Analysts 

To gain some insight into the evolutionary behavior of the 
system, we first consider the case where one of the modes is 
very close to resonance, i.e.. 



(A + e)n. 



(46) 



while all the other modes are non-resonant. We then write 
the orbital and spin frequency evolution equations as 



(47) 





Q tda td 


(Is _ 

fls " 


n! n^ u 1 1 1 

mBa fisfo, tda tsd isi 



(48) 



Here t^ ^ and i^ J are the orbital decay rate and spin up rate 
due to all the non-resonant modesrjThey are approximately 



^ For simplicity, here we do not consider the possibility where 
the star rotates so fast that dynamical tides spindown the star - 
such possibility occurs when the star rotates at a rate somewhat 
beyond Qp (see Lai 1997). 
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given by 



^Ec 



.-1 1 V^o ^ 



m. 



37AB 

TfT T' 

-C/orb 

A£c 



27AJ 



(49) 
(50) 



where A_E = 5I^„ A_Eq and AJ = '^^{rn/aa)/S.Ec, are the 
energy and angular momentum transferred from the orbit 
to the star in the "first" periastron passage (see Section [2|, 
and 7 is the characteristic mode damping rate of the most 
important modes in the energy transfer. Since A J ~ AE/flp 
(where f2p is the orbital frequency at periastron), we find 



fia 



Q' 



O-/ ii L g i Lvj 



(51) 



Note that f^ can be a factor of a few (~ 10) larger than 
^da' ^^^ ^'su is a factor of a few larger than i^^. 

Equations ( |46[ )-(|48| can be combined to yield the evo- 
lution equation for e = Gc/O, — N: 



where 



iV 



1 



mBafls 
NO. 



1 

tsd 



A. 



1 



1 



+ 



td 



(52) 



(53) 



Equation ([52| provides the key for understanding the con- 
dition of achieving mode locking: 



(i) For N > Nc-. The RHS of Equation (521 is always 



negative, and the system will pass through the resonance 
(e = 0) without locking. Physically, the reason is that for 
N > Nc, the orbit decays faster (during resonance) than the 
mode frequency can catch up, so the system sweeps through 
the resonance. 

(ii) For N < Nc and t^^ ^ f^^: Starting from a non- 
resonance initial condition (tin ~ 0.5, or fa ^ 1), the sys- 
tem will evolve toward a resonance-locking state (e = 0) , at 
which 



u 



tda 



> 1 



(54) 



where 5m 
given by 



[Nc/NY - 1. The "equilibrium" value of e is 



1 

2^ 



t^da 



-,1/2 



(7.^)^ 



< 1. 



This resonance-locking state can be achieved when 



(55) 



(56) 



as otherwise resonance "saturation" [/^ ^ (7^^)"^] occurs, 
and the system will sweep through the resonance. Note that 
during the resonance-locking phase, the stellar spin increases 
as 



9.S 



1 

5n 



\fsd ^isu) + 



N^Q. 



NSNinBa^, 
and the orbital frequency increases as 

tl mBa 0.S / -1 -i\ N? 



-t:, 



n nSm ^ 



i^rl Co 



+ 



N^5. 



-t: 



(57) 



(58) 



The above analysis assumes t^^ > 0. In the absence of 



the intrinsic stellar spin-down (i.e 

N, 



^sd = 0), we have 



t~r 



/20fi\ _ 
mN yfi^J ^"^ 



'+t-d\ 



(59) 



Thus i_D can be negative for A^ <, Nc /40. In this case, one 
may expect mode-locking for A^ > Nc- Nevertheless, t^^ is 
only moderately larger (a factor of 10 or so) than f ^^ without 
intrinsic spin-down, so a close resonance with /q ^ 1 (or 
e| <C 1) is unlikely. 



5.3.2 Resonance Locking for m = Modes 

Above (and in other sections of this paper), we considered 
the frequency of a mode in the inertial frame to first order 
in Q,s'. Ga = Ea + mBa^s. In this approximation, m = 
modes have constant frequency and thus cannot lock into 
resonance. However, Lai (1997) find that to second order in 
Q.S , the frequency of m = modes can bo approximated by 






14- 



7 



(60) 



Recall that Sa is the mode frequency in the absence of rota- 
tion. Thus, to second order, the frequency of m, = modes 
can change and resonance locking is possible. 

Let us examine the scenario in which an nia = 2 mode 
with Nm=2 — 90 is locked in resonance, as we suspect is the 
case for KOI-54. Assuming this mode dominates the tidal 
interaction, the system evolves such that 



N„ 



2r2 = (Tq 



(61) 



Then there will be a value of Nm^o for an m = mode such 
that 



N„,=o^ 



(62) 



Using equation (61 1 for the LHS of equation (62l, equa- 
tion ( |60[ ) for the RHS of equation ( |62[ |, and the condition 

o-a,m=o = Nm=o^, we find 



N„ 



12N„ 



l-rUaBa [^/^s — maBa/{2Nni=2)] 



1/2 



(63) 



Using Nm=2 = 90, m^ = 2, Be = 0.84, and fig — ilps — 
16.5S1, we find that Nm.=o — 42. Thus, an tti = mode with 
(ja i; 42n may be able to stay close to resonance for an 
extended period of time. KOI-54 has two very visible modes 
at (T3 — 44f2 and 0-4 = 40n. We speculate that these modes 
may correspond to an m = mode in each star that is nearly 
locked in resonance due to the orbital evolution produced by 
a locked m = 2 mode with Nm=2 — 90. However, because the 
value of Nm^o depends on fig (which is currently unknown 
for KOI-54), we will not consider m = mode locking in the 
remainder of this paper. 



5.3.3 Resonance Locking m Both Stars 

In Section 5.3.1, we considered resonant locking of an m = 2 
mode in the primary star Af . Since the two stars in the 
KOI-54 system are quite similar, resonant locking may be 
achieved in both stars simultaneously. We now consider the 
situation in which an 7n = 2 mode Oa in star M and an 
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m' — 2 mode a^i in star Af ' are both very close to orbital 
resonance, i.e., 



{N + e)n, 



{N' + e')Q, 



(64) 



while all other modes (in both stars) are non-resonant. Here, 
all primed quantities refer to star M' , and unprimed quan- 
tities refer to star M . The orbital evolution equation then 
reads 



n 



tda 



+ 



u 



+ 



id 



(65) 



while the spin evolution is governed by equation ( 57 1 for star 
M and a similar equation for M' . We then find 



e 

N 

./ 
e 



iV 

N' 



-1-x 



where 



fa' /td 
fa/td 



u 

tda 


1 

t^ 


fa' 
tda' 


" t 


Ea' 




Ea ' 





(66) 

(67) 

(68) 



with Ea and E^' the energy dissipation rates (including res- 
onances) due to the resonant modes in star M and A/', re- 
spectively. Thus, analogous to Section 5.3.1, for to > Q and 
t'j^ > 0, the necessary conditions for resonant mode locking 
in both stars are 



A^ < iVc,off, and A^' < K 



c,cff ? 



with 



iVc.eff = 



N, 



y/TTx 



NLa = 



N' 



^r+^ 



37(1 + a;) 



Ba' fJ^a 



1/2 



1/2 



3/'(l-Hx-i) 



(69) 



(70) 



(71) 



Obviously, if the two stars have identical resonances (x = 1), 
then Nc.cS would be smaller than Nc by a factor of \/2- If 
the energy dissipation rates in the two stars differ by at 
least a factor of a few (e.g., x ~ 0.2), then Nc^ch is only 
slightly modified from Nc, while N'^^^g will be a factor of a 
few smaller than A'^^. 



5.4 Numerical Examples of Evolution Toward 
Resonance 

For a given set of modes, the solution of the system of 



equations (42 1- (44 1 depends on the dimensionless param- 



eters tda/tsd, Nc and 7qP, as well as the initial conditions. 
In general, these parameters change as the system evolves. 

5.4-1 A Simple Example 

To illustrate the essential behavior of the secular evolution 
toward resonance locking, we first consider the simple case 
where tda/td and 7q,P are assumed to be constant and iden- 
tical for all modes considered. Figure |4] depicts an exam- 
ple: we include two m = 2 modes (ai and 02), both have 
tda/tsd ~ 10^, ^aP = 10~^ and Ba — 0.8. The parameter 
Nc has an initial value of 110, but we allow Nc to evolve 



via Nc oc f2~^'^. At i = 0, the spin frequency is Q.s — 15f2, 
and the two modes have a\/Q, — 100.5 and a2/f2 = 89.5 
(i.e., both are initially "off-resonance"). We see that the 
stellar spindown first causes cti to lock into resonance at 
ai/fl = 100. The star then spins up, driven by the res- 
onant tidal torque. The second mode passes through the 
A' = 90 resonance, but it produces a negligible effect on 
the ai/Q, = 100 resonance. In the meantime, the orbit de- 
cays, reducing Nc- When Nc becomes less than 100, the res- 
onance locking can no longer be maintained, and ai breaks 
away from the A' = 100 resonance. The stellar spin then de- 
creases (due to the 1/tsd term), which leads to erg to capture 
into the A' = 90 resonance. Eventually, as Nc decreases (due 
the resonant tidal torque) below 90, the second mode breaks 
away from the A' = 90 resonance. This example corroborates 
our analytical result of section 5.3, and demonstrates that 
resonance locking can indeed be achieved and maintained 
for an extended period during the evolution of the binary 
system. 

5.4-2 More Realistic Exam,ples 

Having demonstrated how resonance locking can occur in a 
simplified system consisting of only two oscillation modes, 
we now investigate how resonance locking is likely to occur 
in the actual KOT54 system. We solve the orbital evolution 
equations ( 42 l-(|45[l using the actual values oiea,Qa, and 7a 



found in section|3|to calculate each tda - We include eighteen 
modes in our equations: they are the I = 2, m = 0, m = ±2 
modes for six values of Ea such that the initial frequencies 
of the m = 2 modes range from 40f2 to 170f2. We allow the 
value of Kim to evolve as the values of tJa, SI, and e change 
with time. 

Figures [5] and [6] display the evolution of fi, fig, Nc, and 
the values of aa and Ec, for a selected sample of modes over 
a time span of tens of millions years. We use initial values of 
Qo = 0.95Slobs (here fiobs is the observed orbital frequency 
of KOI-54), Co = 0.84, Qs,o = fips — 16.5f2o, and a spindown 
time of tsd = 3 X 10* years. In Figurelm we have doubled the 
orbital decay rate to account for an equal amount of energy 
dissipation in the companion star. 

Let us start by examining Figure [5] The system quickly 
enters a resonance locking state with the (T3_2 mode, where 
the notation ak.m identifies the mode with the fcth largest 
frequency in our simulation with azimuthal number m. The 
resonance locking lasts for over 100 million years, until the 
o"2,2 mode reaches resonance. The 0-2,2 mode then locks 
in resonance for over 100 million years. During resonance 
locking, the orbital frequency and spin frequency increase 
rapidly. The energies Ea of the resonant locked modes can 
exceed 700 their non-resonant values AEa , corresponding to 
an increase in visibility of over 25 times the visibility during 
a non-resonant state. When no mode is locked in resonance, 
the orbital evolution is relatively slow, with the dominant 
effect being the spindown of the star due to the tsd term. 

For the initial conditions and orbital evolution depicted 
in Figure [5] only two of our 18 modes pass through resonant 
locking phases: two m, — 2 modes with initial frequencies of 
02,2,0 — 103 and 0-3,2,0 — 90. The highest frequency m. — 2 
mode included in our evolution does not lock in resonance 
because it has N > Nc a.t all times. The lowest frequency 
modes do not lock in resonance because they do not satisfy 
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Figure 4. Secular evolution of the stellar spin and binary orbit driven by dynamical tides and "intrinsic" stellar spin-down. Two modes 
are included, both having m = 2, ti^^lt^^ = 100 and 7aP = 0.01. The initial spin frequency is Qs = 150. The top panel shows the orbital 
frequency (solid line) and the spin frequency (dashed line), both in units of their initial values. The second panel shows N^ as defined by 
equation | |38[ l. The bottom two panels show the mode frequencies in units of f2. The two light dotted lines on the third panel also show 
(Ti/Q and 0-2/^. The time (on the horizontal axis) is expressed in units of t^^, the intrinsic spindown time of the star. The system first 
evolves into resonance locking at a\ = lOOO and then a different resonance locking at (T2 = 90O. Note that resonance locking cannot 
occur when aa/f^ > Nc- 



equation (561, i.e., they become saturated before they can 
begin to resonantly lock. This occurs because the value of 
Qa for these high-order modes is small, resulting in a small 
orbital decay rate i^^ for these modes. 

Which mode locks into resonance is dependent on the 
initial conditions. However, as can be seen in Figure |5j the 
resonance with the 0-2,2 niode ends the resonance with the 
as, 2 mode because the 02,2 mode has a larger value of Ea 
and will thus dominate the orbital decay. Therefore, we can 
expect the system to evolve into a state in which the mode 



with the largest value of Ea (and with N < Nc) is locked in 
resonance. 

In Figure [5J the resonance locking with the 0-2,2 mode 
is ended by a resonance with the 04,0 mode. The 04,0 mode 
does not lock into resonance because an m, = mode can- 
not change the stellar spin frequency. However, its resonance 
causes the value of Cl/Q to exceed the value of o-a/o-a for the 
locked mode, so that the system sweeps through the reso- 
nance. In other words, the non-locked resonant mode tem- 
porarily decreases the value of td, decreasing the value of to 
such that resonant saturation occurs for the locked mode. 
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Figure 5. Secular evolution of the KOI-54 system driven by a sample of 18 oscillation modes and intrinsic stellar spindown. Top; orbital 
frequency Q (solid line) and spin frequency ils (dashed line) in units of their initial values. Top middle: the value of A^c (solid line) and 
the mode frequencies 0-2,2/^2 + 15, 0-3^2/^^ + 20, and 0-4.2/^2 + 35 (dotted lines) of the modes depicted in the bottom three panels. Middle: 
the mode energy Ea in units of AEa for the (72,2 mode (solid line) and the 0-4^0 mode (red dashed line). Bottom middle: same as middle 
except the solid line is for the (73 2 mode. Bottom: same as bottom middle except the solid line is for the 0"4^2 mode. The initial spin 
frequency is Us ~ 16. 5f!, and the spindown time is i^d = 3 X 10** years. 
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causing it to pass through resonance. Also note that even 
though the 0-4,0 mode does not lock into resonance, it main- 
tains a relatively large energy [Ea > lOAi^c) for a period 
of over 50 million years. This is indicative of the pseudo- 
resonance locking for m = modes described in Section 
[53:21 

Let us now examine Figure[6] in which the orbital decay 
rate has been doubled (i.e., we multiply Cl by 2) to account 
for an equal tidal response in the companion star. The re- 
sults are significantly different: due to the increased orbital 
decay rate, the initial effective value of Nc,eS has dropped to 
~ 100 (see section 5.3.3) so that the (T2,2 mode can no longer 
lock into resonance. Also, the maximum energy Ea/ /S.Ea of 
the modes is larger in this scenario (because the modes must 
enter deeper into resonance to become locked), so the modes 
can create larger magnitude variations while being locked in 
resonance. 

However, the resonance locking events are shorter be- 
cause the modes are nearly saturated during resonance lock- 
ing, allowing resonances with m = modes or higher fre- 
quency m = 2 modes to easily disrupt the resonance locking 
state. In the bottom three panels of Figure [6] we have plot- 
ted the energy Ea of the 02,0 mode. In Figure [6J three of 
the resonant locking events for the as, 2 mode are ended due 
to a resonance with the (T2,o mode. We have examined the 
results carefully and found that all of the resonance locking 
events were ended due a resonance with an m, = mode or 
an 771 = 2 mode for which N > Nc- Although the resonance 
locking events are brief at the end of the simulation, the 
locking event with the 0-3,2/^2 = 90 mode lasts for about 20 
million years at the beginning of the simulation when the 
system's parameters most closely resemble those of KOI-54. 
We thus conclude it is likely to observe a system such as 
KOI-54 in a resonance locking state. 

We note the secular evolution analysis presented above 
does not take into account the fluctuation of mode ampli- 
tudes due to the changing strength of the tidal potential 
in an elliptical orbit. Therefore, in addition to the results 
displayed above, we have also performed calculations of the 
full dynamic evolution of the stellar oscillations and binary 
orbit, including the back-reaction of the modes on the orbit. 
We find that resonance locking can indeed be maintained for 
extended periods of time, and that non-secular effects have 
little impact on the results discussed above. 



6 OSCILLATIONS AT NON-ORBITAL 
HARMONICS 

As discussed in section 4, in the linear theory, tidally-forced 
stellar oscillations give rise to flux variations at integer mul- 
tiples of the orbital frequency. Many of these oscillations 
have been observed in the KOI-54 system. However, Welsh 
et al. (2011) also reported the detection of nine modes that 
do not have frequencies that are integer multiples of the or- 
bital period, and their observation requires an explanation. 
One possibility is that one (or both) of the stars in 
KOI-54 are 5-Scuti variable stars. The masses, ages, and 
metallicities of the KOI-54 stars put them directly in the 
instability strip and so 5-Scuti-type pulsations are not un- 
expected. However, as pointed out by Welsh et al. (2011), 
the complete absence of any modes detected with periods 



less than 11 hours is troublesome for the (5-Scuti interpreta- 
tion, as most observed (5-Scuti pulsations have periods on the 
order of a few hours. One would then have to explain why 
only long period (P > 11 hours) 5-Scuti-type oscillations 
are visible in this system. 

Another possibility is that some (or all) of the non- 
harmonic modes are excited via nonlinear couplings with 
resonant modes. In particular, it is thought that three-mode 
resonant coupling (parametric resonance) plays an impor- 
tant role in limiting the saturation amplitudes of overstable 
g-modes in ZZ-Ceti stars (Wu & Goldreich 2001) and (5-Scuti 
stars (Dziembowski & Krolikowska 1985). In the KOI-54 sys- 
tem, when a resonantly excited mode reaches sufficiently 
large amplitude, it will couple to two non-resonant, lower- 
frequency daughter modes (see Kumar & Goodman 1996), 
thereby explaining the observed non-harmonic oscillations. 



Parametric resonance can occur when uj^ 



+ ^/3, 



where oj^, is the frequency of the parent mode and uia and 
LOfj are the frequencies of the two daughter modes. The ad- 
ditional requirement that m^ — iria + mp , where m is the 
azimuthal number of the mode, implies that a^ ~ Oa + (Jp. 
Examination of the KOI-54 data given in Welsh et al. (2011) 
reveals that /2 = h + h (25.195 ^Mz = 6.207 ^Hz + 18.988 
/iHz) exactly to the precision of the measurements, where /p 
is the frequency of the magnitude oscillation with the pth 
largest magnitude variation. This provides strong evidence 
that the p = 5 and p = 6 oscillations are due to modes 
excited via parametric resonance and not via direct tidal 
forcing. 

However, no other pair of observed non-harmonic mag- 
nitude oscillations have frequencies which add up to that of 
an observed oscillation. Some of these non-resonant oscilla- 
tions may be due to modes excited via parametric resonance 
that have undetected sister modes. This possibility is espe- 
cially appealing if one considers the scenario in which one 
of the two daughter modes has 777 = (and could thus be 
easily detected for small values of ia), while its sibling has 
777 = 2 (and is thus very difficult to detect for small is). We 
therefore suggest that all non-harmonic flux oscillations are 
produced by nonlinear mode coupling. Deeper observations 
(with better photometry) and analysis may reveal additional 
evidence for nonlinear mode couplings in the KOI-54 system 



7 DISCUSSION 

We have shown that many properties of flux oscillations de- 
tected in the KOT54 binary system by KEPLER can be un- 
derstood using the theory of dynamical tidal excitations of 
g-modes developed in this paper. In particular, our analysis 
and calculation of the resonance mode locking process, which 
is driven by dynamical tides and intrinsic stellar spindown, 
provides a natural explanation for the fact that only those 
modes with frequencies {(Jo) less than about 90-100 times 
the orbital frequency (f2) are observed. Our result suggests 
that the KOI-54 system is currently in a resonance-locking 
state in which one of the stars has a rotation rate such that 
it possesses a m = 2 mode with frequency Oa = 90f2 and 
the other star has a similar mode with a^^i — 91f2 — these 
modes produce the largest flux variations detected in KOI- 
54. Our analysis shows that the system can evolve into and 
stay in such a resonance locking state for relatively long 
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Figure 6. Same as Figure[5] except this evolution doubles the orbital decay rate due to each mode, and the red dashed line corresponds 
to the f72 mode. 
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time periods, and it is reasonable to observe the system in 
such a state. Other less prominent flux variations at lower 
frequencies can be explained by the m = Q tidally forced os- 
cillations, many of which may be enhanced by the resonant 
effect. We have also found evidence of nonlinear three-mode 
coupling from the published KEPLER data and suggested 
that the nonlinear effect may explain the flux variabilities 
at non-harmonic frequencies. 

In our study, we have used approximate quasiadiabatic 
mode damping rates. Obviously, it would be useful to repeat 
our analysis using more realistic mode damping rates, calcu- 
lated with the full non-adiabatic oscillation equations. Our 
calculations of the flux variations of tidally induced oscilla- 
tions have also highlighted the importance of accurate treat- 
ment of non-adiabatic effects in the stellar photosphere. We 
have only considered g-modes (modified by stellar rotation) 
in this paper. Our general theory allows for other rotation- 
driven modes, such as inertial/Rossby modes. It would be 
interesting to study tidal excitations of these modes as well 
as the nonlinear 3-mode coupling effect in the future. 

More detailed modelling of the observed oscillations 
may provide useful constraints on the parameters of the 
KOI-54 system, particularly the stellar rotation rates and 
spin inclinations. On a more general level, KOI-54 may serve 
as a "laboratory" for calibrating theory of dynamical tides, 
which has wide applications in stellar and planetary astro- 
physics (see references in Section 1). 

While we believe that our current theory provides the 
basis for understanding many aspects of the KOI-54 obser- 
vations, some puzzles remain. In our current interpretation, 
the most prominent oscillations (at 90$! and 91il) occur in 
the two diflferent stars, each having a m = 2 mode resonantly 
locked with the orbit. While one can certainly appeal to co- 
incidence, it is intriguing that the strongest observed oscil- 
lations are at the two consecutive harmonics (A^ = 90, 91) 
of the orbital frequency. More importantly, as discussed in 
Section 5.3.3, when both stars are involved in resonance lock- 
ing, the effective Nc (above which resonance locking cannot 
happen) is reduced. For example, using our canonical stel- 
lar parameters (k = 0.05), we find Nc — 131 [see equation 
( 39 1] , but the effective Nc would be reduced to 92 if the two 
stars have identical modes in resonance locking. In the case 
where the resonant energy dissipation rates in the two stars 
differ by a factor of more than a few, the effective Nc of one 
star would remain close to 131, while the effective Nc for the 
other star would likely be reduced to a value below 90. In 
this case, it would be unlikely to find the binary system in 
the resonance-locking state involving both stars with A'^, A'^' 
close to 90, because it is most likely for modes with A'' just 
less than Nc to be locked in resonance (see Section 5.4.2.). 
Thus, we find it likely that the energy dissipation rates are 
nearly equal (to within a factor of 2) in each star in the 
KOI-54 system. 

Finally, we may use our results to speculate on the evo- 
lutionary history (and future) of the KOI-54 system. The 
star with M2 = 2.38M0 and radius R2 — 2.33Rq has an 
age about 3 x 10* years, according to our stellar model 
generated by the MESA code. During the resonance lock- 
ing phase, the orbital eccentricity decreases on a relatively 
short time scale (of order ~ 10* years). The current large 
eccentricity of KOT54 then suggests that resonance locking 
has not operated for a large fraction of the system's history. 



When resonance locking is in effect, the orbital evolution 
time scale is set by primarily by the spindown time scale, 
tsd, and the value of 5m for the resonant mode [see equation 
(58l]. In the future, the stars will continue to expand into 
red giants and the spindown time scale will decrease. This 
will cause the orbital evolution time scale to correspondingly 
decrease, leading to rapid orbital decay and circularization 
of the system. 

While our paper was under review, a paper by Burkart 
et al. (2011) appeared on arXiv. Burkart et al. carried out 
non-adiabatic mode calculations and also attribute most of 
the oscillations to m = modes, although they do not reach 
a definite conclusion on the source of the A' = 90, 91 oscilla- 
tions. They did not consider resonance locking (instead they 
consider the qualitatively different phenomenon of "pseu- 
dosynchronous locks"), and appeared to attribute all reso- 
nances to random chance. They also showed that 3-mode 
coupling (see Section 6) is possible. 
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